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Abstract 

The investigation of contact symmetries of re-parametrization invariant La¬ 
grangians of finite degrees of freedom and quadratic in the velocities is presented. 

The main concern of the paper is those symmetry generators which depend lin¬ 
early in the velocities. A natural extension of the symmetry generator along the 
lapse function N(t ), with the appropriate extension of the dependence in N(t ) 
of the gauge function, is assumed; this action yields new results. The central 
finding is that the integrals of motion are either linear or quadratic in veloci¬ 
ties and are generated, respectively by the conformal Killing vector fields and the 
conformal Killing tensors of the configuration space metric deduced from the ki¬ 
netic part of the Lagrangian (with appropriate conformal factors). The freedom 
of re-parametrization allows one to appropriately scale N(t), so that the potential 
becomes constant; in this case the integrals of motion can be constructed from the 
Killing fields and Killing tensors of the scaled metric. A rather interesting result is 
the non-necessity of the gauge function in Noether’s theorem due to the presence 
of the Hamiltonian constraint. 
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1 Introduction 


The Lie symmetry method is a systematic tool for the study of systems of differen¬ 
tial equations. The importance of Lie symmetries lies in the fact that they can be 
used in order to reduce the order of an ordinary differential equation or to assist 
in the solution of differential systems by reducing the number of independent vari¬ 
ables. The Lie symmetries which leave invariant a variational integral are called 
Noether symmetries. Noether symmetries form a subalgebra of the Lie symme¬ 
tries of the equations following from the variational integral. The characteristic of 
Noether symmetries is that to each symmetry there corresponds a divergence free 
current along with the corresponding charge. The well known conservation laws of 
energy, angular momentum all follow from Noether symmetries. It is well known 
that conservation laws are important tools which can be used for the determination 
of integrable manifolds of a dynamical system in mathematical physics, biology, 
economics and many others 000 00 00 00. 

The determination of the Lie and Noether symmetries of a given differential 
equation consists of two steps, (a) The derivation of the symmetry conditions and 
(b) the solution of these conditions. The first step is formal and it is outlined in 
e.g. mmm- The second step can be a difficult task since the symmetry con¬ 
ditions can be quite involved. Tsamparlis & Paliathanasis in mm proposed a 
geometric method for the solution of the symmetry conditions for regular Hamilto¬ 
nian systems which describe the motion of a particle in a Riemannian space under 
the action of a potential. It was shown that the Lie and Noether symmetries of 
that system are related with the special projective algebra of the underlying space. 
This geometric approach has been extended to the determination of the Lie and 
Noether point symmetries of some families of partial differential equations mama. 

A similar analysis has been established by Christodoulakis, Dimakis & Terzis 
m for constraint Lagrangians quadratic in the velocities, where it was shown that 
the point symmetries of equations of motion are exactly the variational symmetries 
(containing the time reparametrization symmetry) plus the scaling symmetry. The 
variational symmetries was seen to be the simultaneous conformal Killing vector 
fields of both the metric, defined by the kinetic term, and the potential. 

These results have been applied in various areas, i.e. in classical mechanics, 
in general relativity and in cosmology; in each case new dynamical systems which 
admit symmetries were found for instance see [IS [191 [23 [HI 12111221123 [23 [25] 
and reference therein. 

In this work we extend this geometric approach in order to study the Lie 
Backlund variational symmetries (or dynamical Noether symmetries) of constrained 
Lagrangian systems. Specifically, we study the symmetries which arise from con¬ 
tact transformations, the contact Noether symmetries. Contrary to the point 
transformations which are defined in the configuration space manifold, where the 
dynamical equations are defined, the contact transformations are defined in the 
tangent bundle thus they depend also in the velocities [23 EH- Such contact 
Noether symmetries of regular Lagrangian systems with a potential have been 
studied in [28] and it was shown that the corresponding conservation laws follow 
from the Killing tensors of the kinetic metric and the potential. Some important 
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conservation laws of that kind in physics are: the Runge-Lenz vector field of the 
Kepler problem, the Ray-Reid invariant, and the Carter constant in Kerr space- 
time [231 M ED E2]. As it will be shown, in the case of constrained Lagrangian 
systems the results are different from that of [28] in that the Noether contact sym¬ 
metries are related to the Conformal Killing tensors of the underlying space. The 
method we develop can be used and for regular Lagrangian systems. Furthermore 
it is interesting to note that, in the presence of constrains, the gauge function of 
Noether’s theorem does not play a significant role: the presence of the quadratic 
constraint renders it non-essential in finding the conservation laws. The plan of 
the paper is as follows. 

In section [2l we give the basic theory of the Lie Backlund symmetries. In sec¬ 
tion [3j we study the linear contact symmetries of constraint Lagrangian systems 
and we prove the main theorem of this work. The conservation laws which follow 
from the contact symmetries are studied in|4j Furthermore, in section [5] we demon¬ 
strate the main results in two applications, (a) the determination of the variational 
contact symmetries for the geodesic equations of a Biv pp-wave spacetime and (b) 
the derivation of the Ray-Reid invariant of corresponding constraint Hamiltonian 
system. Finally, in section [6] we discuss our results. 


2 Preliminaries 


For the convenience of the reader, in this section we present the basic properties 
and definitions concerning the generalized symmetries. 

Let H = H (x l , u A , u A , be a function which is defined in the space 

Am = |x l , u A , u A , u A j, ...j where x l are n independent variables and u A are m 

dependent variables and u A = ^~r- The function H describes a set of differential 
equations, for n = 1 of ordinary differential equatons and for n > la set of partial 
differential equations. 

We shall say that the function H will be invariant under the action of the 
following infinitesimal transformation 


x 1 = x l + ef (x\ u B ,u B , u^..) (2.1) 

u A = u A +er] A (x',u B ,u*l,u%...) (2.2) 

if there exist a function A [x 1 , u, u^, u^j ...) such as the following condition holds 

m 

[X, H] = XH , mod H = 0. (2.3) 


where X = §§$* + is the generator of the infinitesimal transformation (12.11) . 

m, i- e - 


A = f (x\u B ,u B ,u B j..) di + r/ A (x\u B ,u B ,u B i:j ...) d , 


(2.4) 


and [.,.] is the Lie Bracket. When condition holds, the vector held X is 

called Lie Backlund symmetry of the function H; therefore we conclude that a 
Lie Backlund transformation is a transformation in Am which preserves the set of 
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solutions u A of H (x*, u, u^, in Am- An alternative way to write condition 

(12.311 is 

X [n] H = A H , mod H = 0 (2.5) 

where is the nth prologation of X in the space Am- 

Every differential equation admits as Lie Backlund symmetry the vector field 


Di = di + u.-,d u + u t ijd Ui + ... which is the total derivative operator [IT]. Let X 
be a Lie Backlund symmetry then the vector field 

X = X — f A = (e fc - f k ) 8 k + (V - f k u A k) d uA + - (2-6) 

is also a Lie Backlund symmetry for arbitrary functions /*. Without loss of gen¬ 
erality we could select p = £* and obtain 

X=(r, A -eu A k )d uA . (2.7) 

The generator (I2.7f) is called the canonical form of (12.411 . Since we can always 
absorb the term £ k u k inside the r/ we conclude that X = Z A u B t] ... J 8 u a 


is the generator of a Lie Backlund symmetry. This form of the generator is also 
suitable for the variational symmetries of the action e.g. see |I2] p. 331-333. 

A special class of Lie Backlund symmetries are the contact symmetries defined 
by the requirement that the generator depends only on the first derivatives uj, i.e. 
it has the general canonical form 

X c = Z A (x i ,u B ,u B i )d u A. (2.8) 

2.1 Noether’s Theorem for contact transformations 

Consider now that H = H (f, q A ,q A , q A ) where q A = and H follows from a 
variational principle, that is there exist a Lagrange function L ( t,q A ,q A ) such as 
(L) = 0 where El is the Euler-Lagrange vector. Let X = Z A (t,q B ,q B ) 8a 
be the generator of a contact transformation. Then, if there exist a function 
F = F (t,q B ,q B ) such as [26J 

X [1] L = ^ (2.9) 

dt y ' 

the vector field A is a contact Noether symmetry of the Lagrangian L (t , q A , q A ) [33] 
Where jW is the first prolongation of X, i.e. A'W = X + Z l 8 q i. Further¬ 
more, we have that for every contact Noether symmetry there exist a function 
I = I (: t,q B ,q B ) [MMi 

l{t,q B ,q B ) = Z A (t,q B ,q B ) -F (2.10) 

such as ^ = 0, i.e. the function / (t,q B ,q B 'j is a first integral of the Lagrange 
equations and called a contact Noether Integral. 

In the following section we study the contact symmetries of singular constrained 
Lagrangian systems. 
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3 Contact transformations linear in velocities 

We are seeking variational symmetries of the action 


A 


JdtL = Jdt 


^G^(q)rq u ~NV(q) 


fi,v= 1 ,d 
det G^ / 0, 


generated by contact transformations, that are at most linear in the velocities q fl . 
As it is well known, m describes a singular system consisting of d + 1 degrees of 
freedom. Its dynamics is specified by the Euler - Lagrange equations: 

' H = 2 ^ G ^^ + V = 0 (3-2a) 

S a = q° + T^q u - ^ q a + = 0. (3.2b) 

Owing to the re-parametrization invariance of the action (EH), system EH 
is singular. This fact is best understood by taking the time derivative of (I3.2ap 
and substituting in the result the accelerations q IJ from (I3.2bl) , whence one gets an 
identity (sometimes modulo (|3.2aj) h The final conclusion is that equations (13.211 
cannot be solved for all d + 1 dependent variables, but only for d of them; the 
remaining variable being freely specifiable. If this variable is indeed selected to 
be a specific function of time, one can integrate (in principle) (|3.2bl) to obtain the 
remaining variables, in which case equation (|3.2al) reduces to a relation between 
constants. 

As shown in HZ], in order to reveal all the symmetries that exist by virtue of 
the constraint equation (I3.2ap , N must be considered on an equal footing with the 
configuration space coordinates q^. Thus, it is imperative to avoid any unnecessary 
gauge fixing, prior to the derivation of the symmetries. In this regard, the canonical 
form of generator X of the symmetry transformation under consideration must be 
of the form 

r\ r\ 

X = (£*(*, Q, N) + K K a (t, q, N)q a ) — + Q(t, q, N, q, N) — , (3.3) 

incorporating N as one of the variables. The infinitesimal criterion that the sym¬ 
metry generator must satisfy is condition (12.9p . 

The first prolongation of the generator is 

V' 11 = V + r -X (3.4) 


with 


_ d£ . d . K a , 

“ Tt + Jt {KaQ ] 

= C* + + Co^ + K a, t q a + K% p eqP + K K a 0 Nq a - K K a T® u q^q u + 

^K K a Nq a - N 2 K k ^V^ 
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where ; o = ^ and equations (|3.2b[) have been used, so as to substitute the 
accelerations q a with respect to the velocities q a , N. Note that (El needs not 
a prolongation term in jf- Tl since the Lagrangian is singular and thus does not 
involve the time derivative of the lapse N. 

It is true that 


+ + Jr ** 0 "* 1 + 

+ jrGx,K^q^q X + 1. K Xafi Nq a q X ~ q x (3.5) 

+ -±-K Xa Nq x q a ~ NKfV^ 


(C + K K a q° 


BL 

dq K 


= 2 N^ 


- m*v K + 


- NV K K K rT q a . 


We can easily observe that only up to cubic terms in the velocities appear in the last 
two equations: Nq 2 and q 3 . Hence, it is reasonable to consider adopting a specific 
form for and F. so that no terms exceeding the required order are generated. 
However, this is not a restriction: the existence of any term not complying with 
this selection is trivial, since for any such addition - for example in II - a proper 
term can also be used in F so that they cancel each other out. In this sense we 
define: 


0 = oj(t, q, N ) + (t, q, N)q a + u/ 2) (t, q, N)N (3.7a) 

F = f(t, q, N ) + /« (t, q, N)q a + /$(t, q, N)q a q P (3.7b) 

( 2 ) 

with being symmetric in its indices. As one can see, the gauge function F 

has been assumed to be free of N. The reasoning behind this is that no terms 
involving N should be produced. This acceleration is impossible to be extracted 
from the derivatives of El, since we have already chosen a replacement rule for 
the maximal number of accelerations that we are allowed to use by the equations 
of motion, namely the cp’s. Under the above mentioned assumptions, we get 


BT i j i 

Q dN = ~ - wU - ^G^ifq'Ff - Vu^q a 

i ;( 2 ) 

- ^ G ^rr - v^ 2) n 

i 

=f,t + f,aq a + f, 0 N - fWr%,q»r + -f^Nq a - N 2 f^G^V, 

+ + /$?V + f^Nq a - 2/$T^W + 

- 2 N 2 f^G^v^ + + falo^rr 


(3.8) 


(3-9) 
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where again (13.2bl) have been applied where necessary. 

When relations (E2D, (ESI), EH) and E9D are substituted in the condition 
ehd, one can demand the vanishing of the coefficients appearing in front of terms 
involving different powers of the velocities. This is so because none of the these 
quantities depends on either q^ or N. By working in this way, from the coefficients 
of Nq 2 , we get 

0 + ^K(,a) - = 2f$ + Nf® 0 . (3.10) 

while out of the cubic terms q 3 we derive 

*(.*,> * j = Nf SL>- P- 11 ) 

From the coefficients of Nq it follows that 

fd) _L f(!) = „ 

This equation can be easily integrated to give 

^ = +Ut,q). (3.12) 

Subsequently, the coefficients of the quadratic terms q 2 lead to 

+ K Mt - = JV/g„, + 1 V/®,. (3.13) 


Finally, from the linear terms in q and N, as well as the zero order terms involving 
no velocities we arrive at 


u 

U! 


1 


(2) = --/o 
V ’ 

i, 2N k l f 

-W - -y- K (»a) F - yJ,H 

u> = -y£ K V, K -±f tt + ^fUv>* 


(i) = 

^ NV 


V 


V 


V 


/$ + 


21V 2 

~v~ 



(3.14) 

(3.15) 

(3.16) 


We proceed in order to solve the system. Substitution of (13.1411 in (13.1011 turns 
the latter into 

= ~ Vnv + N fm + l)> (3-17) 

where S, ia is a symmetric matrix. By inserting (13.171) in (13.111) and with the help 
of (I3.15H we deduce that 


V 




which due to ()3.12l) becomes 


V’ a 1 

0) + —^-S<r{nG a p) = 2 
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From this equation we reach the conclusion that - since S a p, G a p, and V are 
not functions of N - the following relations 

i» = Uq) (3.18) 

V' a 

y (3.19) 

must hold. 

At this point, we substitute (13.1611 in (13.13|) and by virtue of (13.1711 we are led 
to 

■c 5 0«„ + + IfVA, - = Ar/g rt (3.20) 

which, by use of (13.1211 . becomes 

If V’ a ~ \ 1 

2 ( G^G a fj, y~^trG a fj, j + — q) = 0. 

Again, due to the fact that none of the remaining functions depends on N, we have 


S af i = S Qll (q ) 

V’ a ~ 

Cf Qil = -—ZoG, 


OLfl ' 


Finally, we are left with the relations 


(3.21) 

(3.22) 


^ = Nf^(t,q,N)+^(q) 

K(«,) = Nf$(t,q,N) - fit 2 £y ] G a , + ^S QfM (q) 


N 


“ = ~ 77 


V 


V 


OJ 


, ( , 1] = --s^V’" - ( L 


V 


V 


“ (2) = -<v 


together with the conditions 


V’ a ~ 

V’ a 

S(otP;ti) = ySa^Gapy 


(3.23a) 

(3.23b) 

(3.23c) 

(3.23d) 

(3.23e) 


(3.24a) 

(3.24b) 


We recognize (|3.24all as the variational Noether (conditional) symmetries found 
in m for singular Lagrangians ensuing from cosmological models. These symme¬ 
tries are generated by vectors f 1 ' 1 in the configuration space, that are simultaneous 
conformal Killing vectors of the mini-super metric G a y and the potential V with 
opposite conformal factors. The second condition, (|3.24bj) . is the one that has to 
be satisfied by a symmetric tensor S a p in order for a contact symmetry generator 
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to exist. As can be seen, S a g must be a conformal Killing tensor with a conformal 
factor that depends on the potential V. 

At this point let us make a comparison between regular systems and re- 
parametrization invariant theories (with quadratic Lagrangians): Regarding the 
former, we know that the relevant conditions require f 11 to be a Killing field of 
both the mini-supermetric and the potential [EJ], i.e. 

Cf ap . = 0 and eV,, = 0 (3.25) 

and, in the case of linear contact symmetries, S a p to be a Killing tensor [28] 
satisfying 

5 \afrvi) = 0 and SjV^ = /, a , (3.26) 

where / is a gauge function. As one can see in each case (either Noether point or 
contact symmetry), two major equations have to be satisfied for regular systems. 
On the other hand, singular Lagrangians require one condition for each type of 
symmetry f (|3.24a|l or (13.24b[) b This allows for an equal or (in most cases) larger 
group of symmetries to be admitted in the case of reparametrization invariant 
systems. A regular and a singular system with the same configuration space metric 
G a /3 and the same potential V, do not necessarily exhibit the same number of 
symmetries; see the relevant discussion in the last section. 

We can reformulate the above conditions (13.2411 in order to show explicitly the 
conformal nature of the them as follows 

= wG a)1 , C^V + LdV = 0 

= SapV’ + ifpV = 0 . 

Additionally, in the case of singular systems one can make use of the re¬ 
parametrization invariance of the theory to simplify conditions (13.2411 . This can be 
done in the constant potential parametrization, where the scaling N i —> N = N V 
is performed. Under this change, the Lagrangian reads 

L=LG^qGr-N, (3.27) 

with Gpv = VGpp being the scaled, by the potential, mini-supermetric. In this 
parametrization, were V = 1, relations (I3.24p reduce to 

Cf a p = 0, S (q/3;m) = 0, (3.28) 

the covariant derivative being now constructed by Gp U and the two tensors S a p,S a p 
related by S a p = V 2 S a g. Thus, the symmetry generators become Killing fields and 
Killing tensors. All the above results can be gathered in order to formulate the 
following theorem: 

Theorem 1. For the constrained Lagrangian in action <ghd and the symmetry 
generator (lH the following statements are equivalent 


9 












1. The vector field is a conformal Killing vector field of the metric G a g and 
the tensor S a p is a conformal Killing tensor of the metric G a g obeying the 
conditions 


CfiG atl = uG ati , C~V + wV = 0 (3.29a) 

S{ a ^) = ^^G a p), S rT f i V’ ,T + if/jV = 0 (3.29b) 

2. The vector field is a Killing vector field of the scaled metric G a p = VG a p 
and the tensor S a p is a Killing tensor of the conformal metric G a p 

Cfaa = 0 , A*fiw)= 0 - ( 3 - 3 °) 


4 Conserved quantities quadratic in veloci¬ 
ties 


We now investigate the conserved quantities that are produced by the generator 
we derived in the previous section. From the previous section we have that the 
general form of the generator X of a contact Noether symmetry for the singular 
Lagrangian m is as follows: 


X=(Nfm-(t,g,N)+e( q ))X 

+ 1,N)~ ifX 



d 

dN 





N 


d 

dN' 


However, as we shall immediately see, the terms involving the arbitrary func¬ 
tions /, and / I'm) are trivial. From ( (12.101) 1 we split the produced quantity / 
into four parts 

I = Q + Iq + 1\ 4- 12 

with 


Q 

h 


h 


h 


c<* + — S a 
* + N b » q ) dq« 


f -a 

q 


2 NV 

NfjpGi** 

NfflG 


dL 

df° 

dL 


dq c 




dL 

dq a 


(4.1) 

(4.2) 

(4.3) 

(4.4) 


It is easy to prove that Iq, I\ and I 2 construct trivial integrals of motion, i.e. 
the corresponding conserved quantities are zero on the constraint surface. 
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• For /(): 


! «-f = - , = - f (jw G " w +') = 0 

due to the constraint equation (|3.2al) . 

• In the case of Ip. 

h ~ = Nf^G^±G aa r - 4V = o. 

• Lastly for 12'. 



Thus, without any loss of generality we can consider the gauge function F to be 
constant and deduce that the only existing conserved quantity is 


Q 


?a dL , oa0 dL dL 

dq a dq a dqP ’ 


(4.5) 


as derived by the only non trivial part of the generator 


X = 




V’V / ~ \ 

— [N^ + 2 Vf) 


d 

8N' 


with £ a and S a /3 satisfying (I3.24al) and 13.24bl) respectively. As seen by (14.51) . in the 
phase space, the vector gives rise to integrals of motion linear in the momenta 
Pa = Jpr, while the tensor S a ^ quadratic. As a result, we can state the following: 

Theorem 2. A singular system described by a Lagrangian of the general form 

L = ^G IM ,{qW<F-NV(q), (4.6) 

admits an integral of motion quadratic in the momenta of the form (p^ = ) 

Q = S a/3 (q)p a pp, (4.7) 

when S a p obeys condition (13.29bl) . 

Conformal Killing tensors (CKTs) can be divided into two major classes: re¬ 
ducible and irreducible. The first type consists of tensors that are trivially con¬ 
structed as tensor products of other CKTs of lower rank. In the case of second rank 
conformal Killing tensors, the reducible CKTs are made up by conformal Killing 
vectors. Of course, if there exists an integral of motion linear in the momenta, 
we expect its square to be again a constant of motion, i.e. if there exist £^’s that 
satisfy (|3.24al) , second rank tensor products can be constructed by them satisfying 
(13.24bD . Indeed, let us consider a second rank tensor S a p = Awith A 1J 
being a symmetric constant matrix and £j, fj vectors that satisfy (I3.24al) (the 
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indexes I, J are just used to discriminate among the various vectors). Under these 
considerations we write 


= -g“ + ^;/i) £/a + + £/>;<*) |j/3 + + fl/3;a) O/, 

+6a (£jp;n + £jn;p) + £ip ( £jor,n + + Zip (' £jar,j3 + £j0;a) 

K IJ V’ (T ~ ~ ~ ~ ~ 

= g [^Ia^jpGaij, + flaCjaG/S^ + flcrfj ^G a p 

J rflaf,J(rGpf 1 + flpfjcrGafj, T flfif,Jo-G a /3 


= ~ A /J — fiafj^G^ 


V’’ 


V ^ 

T 


‘*-kr(/i G C yjj^ 


and thus symmetries defined by £ a, s can be trivially used to construct contact 
symmetries generated with the help of reducible conformal Killing tensors. An 
occurrence quite common in cosmology is that of conformally flat configuration 
metrics in mini-superspace; in this event it is well known that irreducible CKTs 
do exist, see [35J. 


5 Examples 

In order to make things more transparent we apply the previously developed gen¬ 
eral theory to two specific examples; the first is the determination of the geodesics 
of a pp-wave spacetime while the second concerns the constrained two dimensional 
Hamiltonian Ermakov-Ray-Reid system. 

5.1 Determination of Geodesics 

When the potential V(q) in the action (13.11) is constant the corresponding La- 
grangian can be used in order to calculate the geodesics of a manifold with metric 
tensor G MI/ ; equation ()3.2bD represents the geodesic equation in the non-affine pa¬ 
rameter t. Since we have theorem d2j) at hand we can fix the gauge, i.e. choose a 
particular functional form for IV(t); the most natural choice for the time variable 
is N(t) = 1 which gives the time variable t = s the role of an affine parameter. 

In |36] the authors calculated the Killing tensors for a number of pp-wave 
spacetimes along with their conformal Killing vector fields. The metric tensor G 
for the type Biv pp-wave spacetime admits the form 

2 

G = Gapdq 01 <g> d q p = — 2 du <g> dv - du ® dii + dy <g> dy + dz ® dz. (5.1) 

The conformal algebra of the above spacetime is six dimensionsal Sq D D Q\ 
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(four Killing vectors, one homothetic, and one conformal) with basis 


Xi = d v , X 2 = d u , X 3 = d y , X 4 = yd v + ud y 

X 5 = 2 ud u + yd y + zd z , X 6 = u 2 d u + \{y 2 + z 2 )d v + u(yd y + zd z ). 

The five irreducible Killling tensors are 

(s,U = - v%sh 

(SlU = 2yz-H^ - 

(S 3 )^ = 2- IW^J, - 

(SlW = 2 u2 - 2 ^^ ) - zStJh + 

(SsU = (z 2 + 2u 2 z- 2 )S'(Jl ) - 2 uzS^Jh + uH^y 

Since in this case the potential is zero the symmetry generators are the four 
Killing vectors Xi, i = 1,4 and the five Killing tensors Sj, j = 1, 5, which produce 
the nine integrals of motion {/* = ( Xi) a p a , Qj = {Sj) a pp a jr} along with the 


quadratic constraint W (I3.2al) : 

h = —u, h = -2 z~ 2 u — v, I 3 = y, I 4 = yu — yii (5.2a) 

Qi = z~ 2 (~z A y 2 + 2 yz 3 yz - ( 2 u 2 + z 2 z 2 ) y 2 ) (5.2b) 

Q 2 = z~ 2 (2 yu 2 - z 3 yz + yz 2 z 2 ) (5.2c) 

Q 3 = z~ 2 ( z 3 u (zy — zy) + uy (2 u 2 + z 2 z 2 ) — uz 3 yz) (5. 2 d) 

Q 4 = z~ 2 (2 uii 2 — z 3 uz + uz 2 z 2 ) (5.2e) 

Q 5 = z~ 2 (u 2 (2 u 2 + z 2 i 2 ) + z 4 u 2 — 2uz 3 uz ) (5.2f) 

T-L = t ^ 2 (z 2 (2 — 2iiv + y 2 + i 2 ) — 2 ft 2 ) , (5.2g) 


where the dot represents differentiation with respect to the affine parameter s. 
From the integrals I\,I 3 we can solve for u(s),y(s ) i.e. 


u(s) = -h s + c u , y(s) = I 3 s + c y , 


(5.3) 


where c u ,c y are constants of integration. The next step is to eliminate the i 2 from 
the integrals Q 2 , Q 4 using a linear combination of them, yielding 

dz 2 

—-— = J\ + 2J2S => Z 2 = J 2 S 2 + TlS + Cj;, 
as 

where c 2 is another constant of integration and the constants J \, J 2 are combi¬ 
nations of the existing constants, i.e. Ji = 2 ( Q 4 c y — Q 2 c u ) / (c y I 1 + 0 ^/ 3 ), J 2 = 
(Q 2 A + Q4I3) / (cj,/i + c u I 3 ). Finally we are left with the function v(s) which can 
be determinated from the integral I 2 which now assumes the form 


2/i 

J 2 S 2 + Jis + cz 


I 2 =>■ v(s) 


_4 h 

y/ 4 c z J 2 


= arctan 
«/i 2 


2J 2 s + J\ 

V*c z j 2 - J'{ 


I 2 S + Cy, 
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with c v the last constant of integration. If one counts the constants that appear 
in the solution space the total number of them is 9 but the actual number that 
describes the geodesics equation is 2 x 4 = 8, thus there must exist a relation 
between the 9 constants. The desired relation comes from the integral Q 2 and 
reads c z = (8 if + Jf) /(4J 2 ). 

Gathering the above results the final form of the solution space is 


u(s) = -I lS + c u , 
y(s) = I 3 s + c y , 


2 J 2 S + J\ 
2n/27, 
( 2 J 2 S + J\) + 8 if 


v(s ) = \/2arctan 

z 2 ( S ) = 


— I 2 S + Cy 


4 J 2 


(5.4) 


It is interesting to note that, since we have not used the Hamiltonian constraint 


someone insists on the validity of the constraint, then as stated before, a relation 
among the constants emerges, i.e. 


(5.2g) to solve the system, the solution (15.411 describes an unconstrained system; 


<^2+^3 — 2/1/ 2 + 2 — 0. 


5.2 Hamiltonian Ermakov—Ray—Reid system. 

The two dimensional constraint Hamiltonian Ermakov-Ray-Reid system can be 
described by the Hamiltonian 

H = 2n( p °‘ +p p) + NV ’ V = 2 W (“ 2 +^ 2 ) + ^2 J (“) > ( 5 - 5 ) 

with lu, J arbitrary functions of their arguments. The origins of this system can be 
traced back in 1880 in the pioneer work of Ermakov m (see [38] for the English 
translation) along with the generalization of Ray and Reid }39t [31] ; for the uncon¬ 
strained Hamiltonian version see mm- The theoretical interest in this system 
resides in its admittance of an integral of motion for every function u>, J, namely, 
the Ray-Reid invariant 



The above integral can be produced from a Lie point symmetry with the help 
of Noether’s theorem, as first pointed out by Lutzky in |42]. Here we are going to 
show how the integral (15.611 emerges naturally from the contact symmetries of the 
Lagrangian. 

The metric tensor Gy v = Sy U emerging from (15.511 corresponds to a flat two- 
dimensional space. It is well known that the Killing and/or conformal Killing 
tensors span an infinite dimensional space when the dimension of the base manifold 
space equals to two. The CKTs S a p along with the corresponding conformal factors 
ip a are 

a -( F(a,/3) h - | (fi(z) - f 2 {z)) \ , , 

Q/3 \ki-§(h(z)-f 2 (z)) F(a, f3) + fi(z) + f 2 (z)J ’ 1 J 
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and 


V’a = ( d a F(a , /3), d^F(a, /3) + f[(z) + f 2 (z)) , (5.8) 


where z = (5 + ia, z = /3 — ±a and k\ = const. 

In order to apply theorem ([!]) we must demand the validity of the condition 
(13.29bl) where the potential and the CKTs are given by equations (15.511 and (J5.7D re¬ 
spectively. It is interesting that we can satisfy condition (13.29bl) for every function 
u and J when the functions F, fi , f 2 are given by 


F(a,(3) = 


(4c 2 J(/3/a) — 2 c 2 /3 2 ui(a 2 + /3 2 ) + c 3 ) cr 


2 J(f3/a) + a 2 u(a 2 + f3 2 ) 
fi(z) = ci + c 2 z 2 , f 2 (z) = -ci + c 2 z 2 

With the above choice the resulting CKTs are 

(2 a 2 J{w) — a 2 [i 2 (jj{u) a \ 

(*S"l)a/3 = 


2 J(w) + a 2 u:{u) 
a(3 

,,2 


a/3 

2(3 2 J(w) — a 4 u:(u ) 


(S 2 )a/3 = 


a 


2J(w) + a 2 (jj(u) ) 

\ 


2 J(w) + a 2 oj(u ) 
0 


0 




2 J(w) + a 2 uj{u) / 


(5.9) 


(5.10) 


where w = (3/a, u = a 2 + (3 2 . From these CKT, the respective integrals of motion 
along with the Hamiltonian constraint read 


Q i 


Q 2 

n 


2 


(aa + 



2 J(w) + a 2 u(u ) 


^d 2 + /3 2 ^ a 2 

2 J(w) + a 2 oj{u) 

\ (« 2 + + “2 J ( W )' 


(5.11a) 

(5.11b) 

(5.11c) 


Solving the Hamiltonian constraint for the function oj(u) and substituting in Q\ 
we recover the Ray-Reid invariant (ESD. We would like to remark that with a 
similar approach we can construct and the correspoding conservation law of the 
generalized Ermakov system in an n-dimensional Riemannian space [ 145] . 


6 Discussion 

In this work we have generalized previous results in m concerning the deter¬ 
mination of Lie point - Noether symmetries for constraint systems whose action 
is quadratic in the velocities. The generalization concerns the consideration of 
contact symmetries, i.e. of generators which depend linearly on the velocities. 

The key ingredients of the present approach are: 
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(a) the treatment of the lapse function N(t) and q a {t )'s on an equal footing in 
the generator (hi.•Ill , despite the fact that the corresponding Euler-Lagrange 
equation is djyL = 0, which reveals a larger number of symmetries and 

(b) the dependence of the gauge function F(t, q a , q a ) (13.7bl) in the velocities q a 
which in turn makes it constant, i.e. non-essential in contrast to the regular 
case. 


If someone chooses not to implement the IV-dependence in generator (13.311 then 
the function f i(t,q a ,q a ) (I3.7al) should be zero, i.e. ui = 0 ,lo ft = 0,u/ 2 ) = 0. Then 
equation (I3.23df) u = 0 yields / = f(q a ),C^V = 0 and equation (I3.29al) reads 




= 0. Furthermore from equation (13.23dl) = 0 we have the equality 

S^ a V’ a = —V/2 (/ /V) u which forces the conformal factor to be tpy = \ (f/V) u 
due to equation (I3.29b[) . At this stage if we redefine the conformal Killing tensor 
S a /3 via S a p = S a /3 + /V) G^ u then we arrive at S a p = 0. Recapitulating 

the above situation we conclude that if we do not use the IV-dependence the 
symmetry generators are described by the Killing vector fields and Killing tensor 
fields of the metric Gy V and not the conformal ones. Thus if we had at hand a 
metric G^ v that does not admits Killing fields but do admits conformal ones then 
one could conclude that there are no symmetries, which is of course a fault result. 
In order to give an example let us introduce the 2d metric g\ 


g = g a/ 3 dx a <g> dx^ =>■ g = 


x 2 y 2 — 1 


(xdx (8> dx + ydy (g> d y). 


( 6 . 1 ) 


It is an easy task for one to see that metric (16.11) does not admit Killing vectors but 
since it is two dimensional admits an infinite number of conformal Killing fields. 
If we choose the potential V = (y(x 2 y 2 — 1)) then there exist the two conformal 
Killing fields 

6 = -p<9 x , 6 = xd x + ydy (6.2) 

v x 


that satisfy condition (I3.29al) . thus generating two symmetries. Equivalently one 
can state that fixing the gauge might result in loosing symmetries. 

As far as the issue concerning the dependence of the gauge function F in the 
velocities q a , we would like to stress the following facts: Let’s for a moment follow 
the common practice and demand no velocity dependence in F, i.e. F = F(t,q a ), 
repeating the procedure for the evaluation of the symmetries we arrive instead of 
condition (I3.29bl) to the condition 


SayV'° + fiyV + = 0 , 


which might tempt one to think that it is more general, due to the appearance of the 
arbitrary function f(q a )', however, this is not the case: If we redefine the conformal 
Killing tensor S a p in the same spirit as we did before; i.e. S jW = S IW — ^ (f/V) Gy, u 
the above equation reduces to condition (I3.29bp signaling the non-essentiality of 
the gauge function /. As we have shown the function / generates the Hamiltonian 
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constraint Ti, see equation (14.21) . thus one should have expected the non-necessity 
of the gauge function /, due to the presence of the natural ’’gauge function” of the 
problem in hand, i.e. the constraint T~L. 

Furthermore, our results can also be used to the unconstraint case, provided 
that we have enough symmetries to use, so that we do not need the constraint Ji\ 
as, for example, it happens the case of section O where we calculated the geodesics 
of the pp-wave spacetime. 

Finally, we conclude that the determination of the Noether symmetries (point 
and contact alike) of systems with Lagrangian (13.11) . has been reduced to a problem 
of differential geometry; that is, to the determination of the Killing vectors and 
Killing Tensors of conformally related spacetimes. 
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